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Abstract 

The u parameter of the supersymmetric standard model is replaced by AS", where S is a 
singlet chiral superfield, introducing a Peccei-Quinn symmetry into the theory. Dynamics 
at the electroweak scale naturally solves both the strong CP and ji problems as long as 



A is of order y Mz / M p \ or smaller, and yet this theory has the same number of relevant 
parameters as the supersymmetric standard model. The theory will be tested at colliders: 
the u parameter is predicted and there are long-lived superpartners that decay to gravitinos 
or axinos at separated vertices. To avoid too much saxion cold dark matter, a large amount 
of entropy must be produced after the electroweak phase transition. If this is accomplished 
by decays of a massive particle, the reheat temperature should be no more than a GeV, 
strongly constraining baryogenesis. Cold dark matter may be composed of both axions, 
probed by direct detection, and saxions, probed by a soft X-ray background arising from 
decays to 77. There are two known possibilities for avoiding problematic axion domain 
walls: the introduction of new colored fermions or the assumption that the Peccei-Quinn 
symmetry was already broken during inflation. In the first case, in our theory the colored 
particles are expected to be at the weak scale, while in the second case it implies a good 
chance of discovering isocurvature perturbations in the CMB radiation and a relatively low 
Hubble parameter during inflation. 



1 Introduction 



A spontaneously broken global symmetry remains an attractive solution to the strong CP prob- 
lem [1]. The strong CP parameter 9 is canceled by the dynamical relaxation of the resulting 
pseudo-Goldstone boson, the axion [2]. A global symmetry with a QCD anomaly, U(1)pq, can 
be implemented in one of the simplest extensions of the standard model — models with two Higgs 
doublets hip [1]. In order that /„ be much larger than the weak scale, rh, the primary breaking 
of U(1)pq must come from an electroweak singlet scalar, s. As in the DFSZ invisible axion mod- 
els [3], these scalars will have U(l)pQ-preserving interactions such as sh\h 2 or s 2 hih 2 , but not 
s*hih 2 or s* 2 hih 2 . This extension of the standard model fits well with supersymmetry (SUSY), 
except for an immediate question of why the masses of the two Higgs doublets are much smaller 
than the scale f a . This corresponds to the \i problem in the supersymmetric standard model. 

In this paper we point out an extremely simple model that simultaneously solves both the 
strong CP and \i problems without severe fine-tuning. Ours is certainly not the first such theory 
(i.e., SUSY DFSZ axion models) [4, 5, 6], but it is very simple and has important consequences for 
signals at accelerators and in cosmology. In particular, the gravitino and axino (fermionic SUSY 
partner of the axion) are much lighter than the weak scale, so that the superpartners produced at 
hadron colliders end up decaying to either gravitinos or axinos possibly with separated vertices. 
The saxion, the scalar SUSY partner of the axion, is also lighter than the weak scale. All SUSY 
particles around the weak scale are unstable, but cold dark matter may be composed of axions, 1 
and possibly saxions. The cosmological saxions lead to astrophysical X-ray signals. 

The model is presented in section 2. Limits on the SUSY-breaking scale, and hence on the 
gravitino mass, are discussed in section 3. Section 4 is devoted to the thermal history of the 
model, including late-time entropy production required from saxion evolution. The final section 
contains a summary of the predictions of our theory, along with remarks on the difference of our 
theory from those in [4]. 

2 Model 

We consider a supersymmetric theory at the electroweak scale with a superpotential 

W = XSH^z (1) 

1 It is remarkable that the cosmological dynamical relaxation of the axion field during the QCD era can lead 
to the observed amount of dark matter in cold axions, providing the symmetry breaking scale is f a w 10 11 GeV, 
one order of magnitude above the lower bound set by the SN 1987A constraint. This value for f a required for 
the dark matter may be obtained as the geometric mean of the supersymmetry (SUSY) breaking scale and the 
Planck scale [4]. 
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together with Yukawa interactions of the quarks and leptons to the Higgs doublets Hi )2 . This is 
certainly a very simple theory: it is the minimal supersymmetric standard model (MSSM) with 
\i replaced by AS, where S is a gauge singlet superfield, and the soft parameter fxB replaced by 
AXs. Alternatively it can be viewed as the next-to-MSSM without the kS 3 interaction. There is a 
U(1)pq symmetry that is spontaneously broken by (s), leading to a Goldstone boson. Therefore, 
(s) must be larger than about 10 10 GeV to avoid laboratory and astrophysical constraints. One 
immediate problem is to obtain an appropriate scalar potential for s; (s) should be non-zero, 
but should not be infinite. The superpotential (1) alone neither destabilizes nor stabilizes the 
flat direction s. 

This theory has been studied before as a solution to the strong CP problem [7], and a stable 
vacuum with finite (s) was found in the analysis of the scalar potential involving SUSY breaking. 
However, the authors of [7] believed that a fine-tune was necessary to obtain a stable vacuum, 
and that an effective /j, c s = A (s) parameter of order weak scale, rh, is obtained only as a result 
of another accidental cancellation between (s) ^> rh and A<1. The /i problem was not solved. 

On the other hand, it is known [8, 9] that the theory with (1) has a natural solution to the 
\i problem, provided that the soft mass term for s is sufficiently small 

\m 2 s \ < \ 2 m 2 (2) 

for some dynamical reason (see section 3). Once the two Higgs doublets H 12 acquire vevs 
Vi = cos f3 v and v 2 = sin (3 v, the potential for s becomes 

V(s) = -(A\ Vl v 2 s + h.c) + \ 2 {yl + v 2 2 )s*s, (3) 

leading to a stable vacuum with (s) ~ A _1 m and generating an effective \i parameter 

fx eS = As = A cos (3 sin (3 (4) 

of order rh for any value of A, providing a very elegant solution to the \x problem. The crucial 
point is that the scale of \x is set by A, not by (s). Thus, our observation is that the theory 
described by (1) provides a simple simultaneous solution to the strong CP and fi problems. 

Of course, to satisfy astrophysical limits f a ~ (s) ^> rh, the coupling A must be very small. 
We will demonstrate in section 5 that some models can yield A ~ ^Jrh/Mp, so that the axion 
decay constant f a lies around 10 11 GeV, naturally giving rise to axions as cold dark matter. 

The spectrum of this theory consists of states in the chiral multiplet S, in addition to those 
of the MSSM [8, 9]. The saxion mass, m s , is of order Am ~ 100 eV x (A/10 -9 ), and the axino 
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mass is of order 2 ' 3 \ 2 rh ~ 10~ 7 eV(A/10~ 9 ) 2 : these very small masses are a distinctive feature of 
our theory. There is no stable WIMP dark matter candidate at the TeV scale. The axino is the 
LSP, but is so light that it will not contribute to cold dark matter — axions and saxions are the 
candidates for dark matter. 

Various astrophysical processes, that led to the constraints on the axion, also provide phe- 
nomenological limits on the light saxion field. 4 The saxion is light enough to be emitted from the 
interior of horizontal branch stars, and energy loss from saxion emission, which would shorten the 
lifetime of helium-burning stars, sets a limit on the Yukawa coupling of the saxion to the electron. 
The emission is dominated by the bremsstrahlung-like process e~ + 4 He — > 4 He + e~ + saxion, 
and the constraint on the Yukawa coupling is given by [10] 

s(x sm2,3 ) 2 ~ 14x10 " 29 ' (5) 

or equivalently, 

/ a >4sin 2 /3 x 10 10 GeV. (6) 

Saxions are also emitted from SN 1987A, carrying energy away from the supernova. The energy 
loss rate through the saxion turns out to be roughly the same as that through the axion [10], 
so that the astrophysical limit on f a is a little stronger than in conventional axion models. 
These more stringent bounds apparently indicate f a close to 10 11 GeV so that axions necessarily 
contribute a significant fraction of the dark matter, but this conclusion requires further scrutiny 
since our theory requires entropy dilution of saxion field oscillations, as discussed in section 4. 

3 A Light Gravitino 

The stable vacuum of our theory crucially relies on an assumption |m|| < A 2 m 2 . There are 
models of mediation of SUSY breaking with vanishing m|, but radiative corrections to m 2 s 
through the interaction (1) are also of order A 2 m 2 . Thus, one should expect either i) an accidental 
cancellation between tree-level and one-loop contributions, or ii) the SUSY breaking is mediated 
at a low energy scale, so that the one-loop correction is sufficiently small. Let us briefly see how 
low the mediation scale should be. 

2 The light axino is due to a see-saw mechanism: the fcrmion component of the S multiplet has a mass of 
order (A (h}) 2 /(\ (s)) <~ X 2 fh. This can also be understood in terms of symmetries. Since S is the only multiplet 
relevant to the axion, the axino mass must be Majorana. The mass should be proportional to A 2 because the 
supcrpotcntial (1) has a spurious symmetry under which phases of S and A are rotated in the opposite directions. 

3 When the Kahler potential has a non-renormalizable term \S^S\ 2 /M 2 , there is another contribution to the 
axino mass of order (/ /M) 2 to 3 / 2 - For M of order the Planck scale, M p i, it does not exceed 10~ 3 eV, due to the 
upper limits on TO3/2 and f a obtained in section 3 and 4, respectively. Thus, none of the discussion in this article 
is changed. 

4 Production of axinos is sufficiently suppressed as R parity requires that they must be produced in pairs. 
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The tree-level potential of the CP-even scalars shows that the smallest eigenvalue of the 
mass-squared matrix is positive when [8] 



\ 2 (v 2 + v 2 ) 



(7) 

A 



M 2 



We roughly 5 take this limit to be |£| ^$0.2. Since the renormalizat ion-group equation for m| is 
given by 



dm 2 s (n) A 



2 



2— (mi + m 2 2 + A 2 + m 2 s ), (8) 



d log pi 8n 2 
the one-loop contribution to £ is of order 

ei-ioop--^ln(^), (9) 

where SUSY breaking is assumed to be mediated at some energy scale M s . Thus, it follows 
from the vacuum stability condition (7) that to avoid any fine-tuning the "messenger scale" 
Mg is at most one order of magnitude higher than the electroweak scale. Models with such a 
low messenger scale are found in [12]. For larger values of the messenger scale, the amount of 
fine-tuning increases logarithmically. For example, for gauge mediated SUSY breaking models 
with a messenger scale of 10 3 TeV, fine-tuning of order 1/10 is required between £ tree and £i_i 00 p- 

There is no conflict between the requirement of a low messenger scale and the large value for 
(s), which is driven by the soft operators. The soft parameters in the one- loop effective potential 
are renormalized at As, because this is the combination that appears in particle masses. Since 
we are interested in the scalar potential for s with As of order the weak scale, the soft parameters 
are evaluated at the weak scale, where they are local no matter how low the messenger scale is. 
The scalar potential of the s field is essentially given by physics at the electroweak scale, even 
though the vev (s) f a is much larger than the electroweak scale. 

The fine-tuning argument above favors a low mediation scale, but does not directly constrain 
the fundamental scale of local supersymmetry breaking v^susy- However, in supergravity the- 
ories all scalars fields typically 6 acquire a supersymmetry breaking mass, giving a contribution 
to m| of order [m 3/2 = F SVSY /V3M pl } 2 , where M pl ~ 2.4 x 10 18 GeV is the Planck scale. For 

5 Detailed analysis using tan/3 > 2.5 and fj, > 120 GeV leads to —0.15 < £ < 0.12. For the effects of 1-loop 
corrections to the scalar potential, see [11, 7]. 

6 When the Kahler potential has a certain form, the ordinary gravity-mediated supcrsymmetry-breaking masses 
are absent, and m s acquires only an anomaly-mediated piece, which is of order (A 2 /47r)aim^ 2 . In this case, 
although to 3 /2 of order the weak scale is allowed, the discussion in the following sections is not modified essentially. 
In section 4, the gravitino is no longer light, but the entropy production required to dilute the saxion oscillation 
also dilutes the gravitino number density, and there is no gravitino problem. It is known that the Affleck-Dine 
mechanism works for baryogenesis [13]. In section 5, we still expect separated vertices in colliders, although they 
arise from decays to axinos rather than to gravitinos. 
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this contribution to satisfy (2) without any fine-tuning, the bound on the scale of local super- 
symmetry breaking is 



If this bound is saturated we would normally expect a gravitino problem: the gravitinos 
are in thermal equilibrium at the weak scale, and although they are somewhat diluted by later 
annihilations, they still give too much hot dark matter. This would lead to an even stronger 
bound on V Fsusy than given above. Although the gravitino is not the LSP, it decays to axion- 
axino with a lifetime longer than the age of the universe, so that the gravitino problem is not 
alleviated. However, in section 4 we see that entropy production after the electroweak phase 
transition is required to dilute saxion oscillations, and this will also dilute the gravitinos. 

4 Thermal History 

4.1 Saxions and Late-Time Entropy Production 

Supersymmetric axion models always involve a saxion field with a mass at most of order the 
SUSY breaking scale, m. Thus, there is a flat direction, and its evolution in the early universe 
must be examined carefully. 

During inflation the saxion field could be zero or large, for example of order the Planck scale, 
depending on its coupling to the inflaton. We begin by supposing that it is at the origin. In this 
case it stays at the origin until the Peccei-Quinn phase transition is triggered by the Higgs vev at 
a temperature of order the electroweak scale, as seen from the potential of eq. (3). Immediately 
after the phase transition, the saxion field oscillates about the minimum of its potential, with an 
energy density, V s , of order m 4 . If the saxion were to decay rapidly enough, for instance, with 
a decay rate of order T ~ m 3 /f£, the field energy would rapidly convert into radiation giving 
no problem. However, the saxion mass is not of order m, but Am. The lifetime of the saxion is 
of order r ~ 10 2 (10~ 9 /A) 5 x 10 10 yrs., and is much longer than the present age of the universe. 
The oscillation of the saxion field, which behaves like matter, overdoses the universe. To avoid 
this we study the dilution of the saxion field oscillations by large entropy production after the 
electroweak and Peccei-Quinn phase transitions. 

Let us suppose, for simplicity, that the entropy is produced via the decays of a massive 
particle, X, which could be the inflaton, curvaton or flaton. During and after the electroweak 
phase transition the universe is dominated by X, and is therefore matter-dominated. While the 
Hubble parameter is much larger than the decay rate of the X particle, the energy density of 



m 3/2 



100 eV x (A/10" 9 ) 




(10) 
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X, p x , scales as oc 1/a 3 , where a is the scale factor. Some X particles decay at a time much 
less than the X lifetime, producing entropy. The total energy density of radiation p 7 ~ T 4 does 
not scale as oc 1/a 4 , but rather as oc l/a 3//2 , because of the continuous entropy supply from 
the X-particle decays. These X decays to radiation clearly dilute the saxion field oscillation 
energy density. This dilution continues until the age of the universe becomes comparable to the 
lifetime of the X particle, when p 7 is also comparable to px (see e.g. [14]). A long X lifetime, 
and therefore a low value for the reheating temperature Tr, leads to more dilution. Any initial 
thermal saxions are diluted to a negligible level, while the current number density of cold saxions 
in the saxion field oscillation is given by 

n, ( T R \ 5 /l00GeV\"/100eV\ 

where T PQ is the temperature of the Peccei-Quinn phase transition, and the dimensionless 
coefficient c is given by 

Here, the last factor U s /p 7 is evaluated at the epoch of the Peccei-Quinn phase transition, and 
9s(Tpq) and gs(T R ) are the effective statistical degrees of freedom when the temperature is 
around T P q and Tr, respectively. The saxion field oscillations contribute to the present energy 
density an amount 

n.fc. = x 10-' x (_^) ( ] , (13) 



where 



c-~(i+?iU?V 



22/ \5J 



n ^m§ * 10 , 5 

7 2.73K 



9s(T R ) \ 



The saxion energy density does not depend on the choice of f a . ^cdm^ 2 ~ 0.1 requires 
Tr ^ 1 GeV, and when this bound is saturated, the saxion is also a significant component of the 
CDM. 7 

Such a large entropy production not only dilutes saxions to an acceptable level, but also 
dilutes other species. The SUSY particles of the MSSM sector annihilate quickly to lighter SUSY 



7 To be more precise, the vev's of the Higgs fields are determined by a thermal potential, and they change 
as the temperature falls. The minimum of the s field is also changing accordingly. Thus, even if the entropy 
production dilutes the energy of the saxion oscillation when most of the X particles have decayed, one has to 
further make sure that such readjustments after the entropy production do not release too much oscillation energy 
for the CDM. It turns out that the Higgs field values are close enough to the vacuum values when the temperature 
is around 1 GeV or lower, and the saxion oscillation due to this late-time readjustments is not cosmologically 
important. 
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particles, which eventually decay to the gravitino or axino. Axino number-changing reactions 
have already decoupled by the electroweak scale, and gravitino number-changing reactions freeze- 
out well before T R , so that both axinos and gravitinos are significantly diluted by the entropy 
production. Hence, even gravitino masses that saturate the bound of (10) do not lead to an 
amount of hot dark matter in conflict with observation. Furthermore, gravitinos and axinos 
provide negligible contributions to the effective number of neutrino generations during BBN and 
CMB eras. 

Entropy production from X decays also dilutes the baryon asymmetry, by a factor (Tr/Tb) 5 , 
where Tb is the temperature at which the baryon asymmetry is generated. This severe dilution 
implies that the baryon asymmetry cannot be created above the weak scale. One possibility 
is that a baryon asymmetry of order unity is created at the weak scale. Another is that the 
observed small asymmetry is created in the out-of-equilibrium decays of the X particles; see, 
[15] as an example, where .R-parity- violating interactions W = UDD are required, and T R has 
to be less than 1 GeV. 

4.2 Relic Axion Energy Density 

Let us now turn our attention to the energy density carried by the axion field. There are two 
significant components: one from axions emitted by axionic strings, and the other from the 
misalignment of the initial value of the axion field from that of the true potential minimum. 
The energy density from the axion phase relaxation is known to be [18] 

0,0 Xl 0-(^f 7 (^p (15) 

in the absence of the entropy production from X decays. This calculation assumes a radiation- 
dominated background at the epoch of relaxation, when the effective mass of the axion m a eff ~ 
0.1m a x (Aqcd/T^) 3 ' 7 becomes comparable to the Hubble parameter. If T R is around its upper 
bound ~ 1 GeV, then this assumption is justified. Given the lower bound on f a at the end of 
section 2, which is more stringent than in conventional models, in our theory the axion energy 
density accounts for a significant fraction of CDM, and axions are clearly a natural candidate for 
CDM. When the decay temperature is well below 1 GeV, massive X particles affect the axion 
energy density in two ways [19] . The extra energy density from X particles delays the epoch of 
the relaxation of the axion field, and the entropy from X decays dilutes the axion energy density. 
Combining both effects, the axion energy density is given by 

^'"^iimHiivf ■ (16) 
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and in order to have axion CDM we require that f a be raised to 



1 GeV 



) 



1.3 



f a ~ 10 11 GeV 



(17) 



Imposing T R ;C 1 MeV from big bang nucleosynthesis, the upper bound on f a is f a £ 10 15 GeV 



If s is large during inflation then the axionic strings from the PQ phase transition are inflated 
away, and the axion energy density is from the misalignment of the initial phase. However, if s 
vanishes during inflation, axionic strings are formed after the PQ phase transition is triggered at 
the electroweak scale. After further cooling during the X dominated era, the energy density in 
strings reaches a fraction (f a /M p \) 2 \ii(f a /H) of the total energy density. When the temperature 
falls so that m afi s becomes comparable to the Hubble parameter, axionic domain walls emerge 
and the string network turns into the boundary of the domain walls. This string/domain wall 
system rapidly disappears by radiating axions, 8 and the number density of axions is fixed: 



where A ranges from \a(f a /H) to unity. This large uncertainty in A corresponds to the disagree- 
ment between [16] and [17] about the typical energy of axions emitted from the string network. 
Since the axion number density from misalignment is also Hf% at that epoch, the resulting relic 
density from strings is A times that of the misalignment axions, just as in the case of radiation 
dominance. Note also that another uncertainty in the relic density arises from assuming that 
the energy density of the string network is converted into axion particles when m a ^& is equal to 
the Hubble parameter. 

In the case that s vanishes after inflation, we have assumed that there is only a single vacuum 
in the direction of the axion field. However, the model presented in section 2 has three vacua, 
due to three families contributing to the U(1)pq[SU(3)c] 2 anomaly. The number of vacua can be 
reduced to one by introducing extra colored particles. For instance, introducing two vector-like 
pairs of chiral multiplets $j(3) + 3>£(3*) (i = 1,2), with a coupling 



gives an anomaly coefficient of 1, so that there is a unique vacuum. If the couplings Aj (i — 1,2) 
are much smaller than A, then these extra particles have already been excluded by data. On the 
other hand, if Aj are larger than A, the SUSY-breaking masses 9 of the extra particles contribute 

8 There must be only one vacuum in the phase direction of s. Otherwise, such string/domain wall system 
cannot disappear. 

9 This argument does not apply if A, are so large that their masses due to (19) are larger than the "messenger 
scale" M s . 



[19]. 



^a,string ~ Hf^A, 



(18) 



AW = XiSQ&S, 



(19) 
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to £i_i 00 p a t order A 2 /A 2 , requiring excessive fine-tuning. Thus, A, ~ A, and, as a consequence, 
these extra colored chiral multiplets are expected at the electroweak scale. If these vector-like 
particles have the same electric and PQ charges as the up or down-type quarks, then they can 
decay by mixing with the known quarks. 

4.3 Early Breaking of PQ Symmetry 

It has been assumed so far in this section that the s field vanishes after inflation until later 
times. If s takes a very large value s during and after inflation, the radial direction of the s 
field, i.e., the saxion field, starts to oscillate when the Hubble parameter is comparable to the 
curvature of its scalar potential. The energy of the saxion oscillation has to be diluted in this 
well, by an amount that depends on s . If sq is of order of the Planck scale, even the 
entropy production from the X decays is not enough. Indeed, p s /px ~ (so/M p \) 2 when the s 
field starts to oscillate, 10 and for X decays to provide sufficient dilution for saxions, we find that 
the following condition should be satisfied: 



SO < SO.max = 10~ 3 M pl X ^ (20) 

The initial phase of s is well-defined and almost homogeneous inside the present horizon 
provided the Hubble parameter during the inflation Hj is smaller than so; the phase fluctuation 
that might be generated during inflation is of order 59 ~ H I /(27rs ) < 1. The initial phase is 
preserved in classical evolution of the s field until the QCD phase transition, when the potential 
in the phase direction emerges. Since the whole universe inside the horizon falls into the single 
vacuum, there is no problem of domain walls even in the absence of the extra colored particles. 
The estimate (15) or (16) for Q m i S . applies to this case, except that i) the initial phase of sq 
is used instead of the average of random phases and ii) the normalization of f a can be 

different because of the different number of vacua in the phase direction [18]. Note that axionic 
string network is not formed after inflation in this case, and the misalignment of the initial phase 
is the only source of the cosmological axions. 

The phase fluctuation 86 leads to isocurvature density perturbations in the axions and radi- 
ation [20], giving 

UAsocurv 27rM pl VMeV , 1 ] 

Since the current CMB data is consistent with purely adiabatic density perturbation, this implies 
H T < 10 11 GeV,/ ^ 30 , which is quite a non-trivial constraint on many models of inflation. 

V R •^max 

If this bound is saturated, the isocurvature perturbation may be observed in future CMB data. 

10 Here, the gravity-mediated quadratic and/or gauge- mediated logarithmic SUSY-brcaking potential for s is 
assumed. 
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5 Predictions and Conclusions 



In this paper we have made the observation that by promoting the \i parameter of the minimal 
super symmetric standard model to a field, both the strong CP and /x problems are solved. How 
will we know whether this theory is correct? 

The first, and most important, test that the theory must pass is that the three parameters 
fi, A and tan f3, which are independent in the MSSM, must satisfy the relation of equation (4) 
[7] . This signals that the electroweak symmetry breaking sector of the theory is governed by the 
superpotential interaction \SHiH 2 , without the S* 3 interaction of the next-to-MSSM, and that 
the vacuum is the one with (s) = fn/X that occurs when \m 2 s \ < \ 2 v 2 [8, 9]. While verification of 
this relation will show how the /i problem is solved, it is clearly insufficient to demonstrate that 
there is a PQ solution to the strong CP problem. For example, it could be that A « 1CT 3 and 
that small explicit symmetry breaking terms give the would-be axion a mass of order a GeV [9], 
so that the strong CP problem is not solved. On the other hand, observing a very small value for 
A would provide a strong indication that (s) = fh/X is large and is consistent with astrophysical 
constraints on the PQ solution. Such evidence for small A could be found in the cascade decays 
of the superpartners produced at hadron colliders, as we now discuss. 

Recall that all the superpartners have masses of order m, except for the axino and the 
gravitino which are much lighter. An important question is the decay mode and decay rate of 
the lightest superpartner amongst those that have masses of order m, the LSP'. For large A the 
dominant decay of the LSP' will be to the axino, a. For example, if the LSP' is a neutralino, it 
would decay with a rate A 2 m to either ha or Za, leading to the spectacular events discussed in [9]. 
The last decay process of the cascade chain takes place well inside the beam pipe. However, in 
axionic theories we have A ^ 10~ 9 and a low value for V^susy so that LSP' decays to gravitinos, 
G, have a branching ratio comparable to or higher than that of decays to axinos. The decay 
rate to the gravitino is m 5 /F^ VSY , and hence 



When V -Fsusy <C 300 TeV, the LSP' decays are mainly to gravitinos. The decay vertices are 
still within the beam pipe, and the event configuration of decays to gravitinos is quite similar 
to that of decays to axinos. Thus, it may be difficult to distinguish this axionic theory with 
V -Fsusy <^ 300 TeV from the theory with A ~ 10~ 3 in [8, 9]. In cases with photino LSP', however, 
Br(LSP' -> a + 7 )/Br(LSP' -> a + Z) is different from Br(LSP' -> G + 7 )/Br(LSP' -> G + Z). 
Determination of the mixing in the neutralino sector might be able to exclude one of the two 
possibilities above. When v^susy is close to its upper bound of 300 TeV, the LSP' decay to 




(22) 
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a gravitino occurs at vertices separated from the primary vertex. 11 While this is also a typical 
signal of gauge- mediated SUSY breaking with such a value of -F S usy, decays at separated vertices, 
combined with a positive test of the relation (4) for /x, would give a strong indication that A is 
small, and that our proposal for simultaneous solutions of strong CP and /x problems is correct. 

It may also happen that the QCD-charged extra particles introduced in section 4 are within 
the reach of hadron colliders. 

If our theory is correct, hadron colliders will tell us that the LSP is not in the MSSM sector. 
They will also tell us, from the LSP' decay rate, that the gravitino mass is too small for gravitino 
cold dark matter. The axion will be the natural remaining candidate for cold dark matter. 

When f a saturates the lower bound 10 11 GeV, and Tr ~ 1 GeV, the cold dark matter consists 
of both saxions and axions. The standard axion dark matter search [21] can detect the axion 
of this theory. The saxion can be detected indirectly from the soft X rays produced by its 
subdominant decay mode s — > 77. There are two sources of this X-ray flux. One is from the 
saxions distributed uniformly throughout the universe, and the other arises from saxions which 
have fallen into the gravitational potential of clusters of galaxies (c.f. [22]). The former is 
observed as an isotropic flux with a continuous spectrum, because the X rays emitted long ago 
and far away appear red-shifted. A rough estimate for the photon flux is given by 



for E 1 < m s /2, where matter dominance is assumed in the first line, and the number density of 
the saxions in (11) is used in the second line. This is the flux predicted in extragalactic space 
- the flux on Earth is reduced by absorption in the Galaxy, and will have a modified energy 
spectrum. Using ROSAT data, the observed extragalactic soft X-ray background is found to be 
30-65 keV/cm 2 /sec/str/keV for the 1/4 keV energy region [23]. Given the uncertainties in the 
predicted flux of (23), the X-ray flux from saxion decays is consistent with observation. We note 
here that the saxion mass m s could be as high as 1 keV in a parameter region, although we have 
typically considered m s ks 100 eV. 

Observations with high angular resolution suggest that 60% (and perhaps more) of the ex- 
tragalactic X-rays can be attributed to the flux from discrete sources, such as AGN's, in the 

11 In the previous section we pointed out that one way to accomplish baryogenesis at low temperatures is by 
introducing large R parity violating operators of the form UDD. In this case the LSP' is expected to decay 
dominantly via these interactions, so that the separated vertex signal is replaced by the signals of R parity 
violation. 




(23) 
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1-2 keV energy region [24]. But it is not clear whether all of the soft X-ray flux is accounted 
for by such sources in the sub-keV energy region [25, 26, 24, 23], and there is still room for 
extra fluxes with particle-physics origins. Further observation with high angular resolution and 
long exposure time will certainly help determine the purely isotropic extragalactic component 
while removing foregound contamination and contributions from discrete sources. Furthermore, 
our X-ray signal from CDM saxions has a different spectrum from those of the discrete sources 
identified in [24]. Those fluxes decrease with increasing energy and have tails that extend above 
1 keV [24], while our signal's flux increases with energy until it is sharply cut off at an energy of 
half the saxion mass. Thus, observations with high energy resolution [27] will help identify the 
X-ray signal from CDM saxions, when combined with a better understanding of the foreground 
absorption and emission. 

The X-rays from saxions bound to a galactic cluster produce a "line spectrum" , with a width 
AE^/E^ ~ 0(v/c) ~ 1CT 2 , and an energy m s /(2(l + z)), with z the redshift of the cluster. The 
photon flux at the peak energy is roughly given by 



G^rlnstnr 1 1 



dVt dE 1 4vr 10- 2 £ 7 



/ ^IclusterfXs -> 77), 



2 



si 3x 10 /cm /sec/str/keV 

/ VL s h 2 \ ( cluster size A / n B | c i ust er \ / m s c 2 \ ( 10 11 GeV \ ' 
X \5Q B h 2 )\ IMpc j v 10- 3 /cm 3 ; ^lOOeVy ^ f a ) 

where the integration in the first line is along the line of sight. In order to identify these line- 
spectrum X rays, observations with both high angular resolution and high energy resolution are 
required. 

As Tr is lowered below 1 GeV, f a increases above 10 11 GeV and these axion and saxion signals 
become too difficult to see: the axion detection rate is proportional to A ~ 1/ (m a f 2 ) [21], and the 
extragalactic photon flux from saxion decay oc A 6 ' 8 , where eq. (17) and E 1 oc m s as Av are used. 
However, the lowest value of f a (and hence the largest of A ~ 10~ 9 ) is the most theoretically 
well-motivated. So far in this paper, the extremely small value for A has been put by hand, but 



it can be obtained naturally as ym/M p i ~ 10 9 , where M p \ is the Planck scale. Indeed, one can 
think of a theory with a superpotential 

W = ^S'S Hl H 2 + ±S'\ (25) 

where M ~ M' ~ M p \. With a negative SUSY-breaking mass-squared of order — m 2 for S', one 



can easily see that (S') /M ~ yrh/M p \. The mixing between S and S' is so small in this theory 
that the phenomenological analysis given in this article with an effective coupling A = (S') /M 
is completely valid. 
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The theory with the superpotential (1) is quite similar to those [4] with 

W = ^S'SH^ + ^S n S'^ (n + 0, 4) (26) 

at first sight, but these two classes of theories are quite different. In our theory (n = 0), S' 
is neutral under the Peccei-Quinn symmetry, and the mixing between S' and S is quite small. 
Thus, the chiral multiplet S is virtually the only one responsible for the spontaneous Peccei- 
Quinn symmetry breaking. This is one of the most important reasons why the axino, which is 
the LSP, is extremely light in our theory (c./. [28]). Another important aspect of our theory 
is that the stabilization of s results only after the electroweak phase transition, so that the \i 
parameter is predicted in terms of tan/5 and the A parameter. Therefore, our theory is not 
merely a particular case of the theories in [4], but is essentially different. Indeed, the theory of 
the invisible axion presented in this article has several predictions that can be tested in the near 
future. 
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